Application of Relativistic Coupled-cluster Theory to Heavy Atomic Systems with 
Strongly Interacting Configurations: Hyperfine Interactions in 207 Pb + 
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This work presents a first time accurate calculation of the magnetic dipole hyperfine structure 
constants for the ground state and some low-lying excited states of Pb + . By comparing different 
levels of approximation with experimental results, we demonstrate the importance of correlation 
effects which reach beyond lower order relativistic many body perturbation theory. Employing 
relativistic coupled-cluster theory we obtain a quantitative understanding of the core-polarization 
and correlation effects inherent in this system and observe completely different trends compared to 
Ba+. 



Coupled-cluster theory has been used to study a 
wide range of many-body systems^]. Although the 
non-relativistic version of this theory has been very 
successfully applied to a variety of light atoms and 
molecules Q, its extension to the relativistic regime 
is rather recent 0, There have been relatively few 
theoretical studies of properties of heavy atomic systems 
based on the relativistic coupled-cluster (RCC) theory. 
Pb + (Z=82) is the heaviest atomic ion that has been 
trapped and cooled so far0, |(|. The magnetic dipole 
hyperfine constants have been measured for the &p 2 P\/2 
and 6p 2 P 3 / 2 states of this ionQ and these data can 
be compared with calculations of the corresponding 
quantities using RCC theory. Such comparisons would 
indeed constitute an important test of this theory. The 
non-linear RCC in the singles and doubles approxima- 
tion with partial triples added in some cases has yielded 
results to an accuracy of about one percent for atoms 
and ions with a single s valence electron 0, U ITo| . 
However, the correlation effects in Pb + are expected to 
be much stronger as it has a 6p valence electron and two 
6s electrons in its outermost core orbital. 

The hyperfine structure constant (4)for the atomic state 
| JM) can be expressed in terms of a reduced expectation 
value 



with /Lijv being the nuclear magnetic moment and 
the Lande's nuclear g-factor (gi). T^ 1 ' can be written as 
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where rj is the radial position of the j electron, ctj is 
the Dirac matrix and Y^q is a vector spherical harmonic. 

We have used the RCC theory in to obtain the atomic 
wavef unctions. As pointed out in our earlier work [T^ 
coupled-cluster theory is equivalent to all order many- 
body perturbation th eory (MBPT). In the open-shell 
coupled-cluster theory [p| [HI the many-body wavefunc- 
tion for a system with single valence electron can be writ- 
ten as 
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where is the creation operator corresponding to a 
valence orbital V and | <E»o ) is a closed-shell determinan- 
tal state built from occupied Dirac-Fock (DF) orbitals. 
T- and S v - are the closed and open shell excitation 
operators respectively. In this work both T- and S v - 
operators are truncated beyond double excitations and 
triple excitations are added on the leading order MBPT 
level. 



Explicitly, the T- operator is defined as 
T = Ti+T 2 
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P3/2 


d 3/2 


d5/2 


fs/2 


fr/2 


g7/2 


g9/2 




38 


35 


35 


30 


30 


25 


25 


20 


20 


Active holes 
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4 
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3 


3 


1 


1 








Active particles 


7 


9 


9 


8 


9 


7 


7 


7 


7 


Upper energy limit (a.u.) 


2800 


2850 


2850 


510 


510 


22.6 


22.6 


22.6 


22.6 



TABLE I: Description of total number of basis functions, active holes and active particles involved in this calculation 



where t v a and t pq b are the amplitudes of the single and 
double excitations from the closed-shell core. Similarly, 
the open-shell excitation operator (S v ) is defined as 



S\ v + S2v 



ala v s p v + 2 X! a t a l a aa v s p , q a (0.5) 

p^v a,pq 



with s*J and being the single and double excitation 
amplitudes involving the valence electron. 

In coupled-cluster theory the expectation value of any 
operator can be expressed as |9j 
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The above expression was applied to compute the 
hyperfine structure constant A ' as given in eqn. l|U.lfl . 

The orbitals used in the present work were constructed 
as linear combinations of Gaussian type orbitals(GTOs) 
of the form ^(| 



k „-c 
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where k = 0,1,... for s,p,.. type orbital symmetries re- 
spectively. For the exponents, the even tempering condi- 
tion 



o-, 



a P l 



(0.8) 



was used. The occupied orbitals are the DF single par- 
ticle states for closed-shell Pb ++ . The virtual V 



N-l 



or- 



bitals [171 [18j were constructed from the closed-shell po- 
tential of Pb ++ using the same Fock operator. All or- 
bitals were generated on a grid using a two-parameter 
Fermi nuclear distribution approximation given by 



P = 



Pi) 



1 + e (r-c)/o 



(0.9) 



where the parameter 'c' is the half-charge radius, and 'a' 
is related to the skin thickness which is defined as the 
interval of the nuclear thickness which the nuclear charge 
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FIG. 1: Goldstone diagrams representing core-polarization in 
MBPT 



States DF approximation Experiment [6] 



6P1/2 


11513.5 


13000 


6P3/2 


918.3 


583(21) 


7S1/2 


7822.9 




Tpi/2 


1983.1 





TABLE II: Dirac-Fock and Experimental results for magnetic 
dipole hyperfine structure constants of Pb + in MHz 



density falls from near one to near zero The number of 
basis functions used for generating the occupied and 
virtual orbitals are given in table I. The active virtual 
('particle') and core ('hole') electrons considered for 
the coupled-cluster calculations are also displayed. The 
upper energy limits above which the virtual orbitals 
were truncated during the RCC computation are given 
in atomic units. We have chosen ao as 0.00825 and (3 
as 2.73 for all the symmetries for the generation of the 
GTO's. 

In table II we present the results for the hyperfine con- 
stants using the DF approximation and compare them 
with the experimental values for the Qpi/2 and 6P3/2 
states - the only two states on which measurements 
have been made. The poor agreement of the results 
indicate the importance of correlation effects for these 
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Lowest order 


6pi/ 2 


7pi/2 


6p 3 /2 


TSl/2 


MBPT terms 


state 


state 


state 


state 


U(D^ J 


11513.5 


1983.1 


918.3 


7000 r\ 

7822.9 


RMBPT(2) 


15722.8 


2578.4 


302.9 


12663.9 


core-pol. 


1506.2 


82.1 


-814.6 


1624.21 


pair-corr. 


2297.4 


359.6 


203.6 


3012.7 



TABLE III: Second order MBPT results for Pb+ hyperfine 
structure constants in MHz (second row) and the dominating 
contributions (third and fourth row) as shown in fig.l D 

Virtual orbital Core orbital RMBPT(2) RCCT 



7si /2 

8S1/2 
9Si/2 
10Sl/2 
HSl/2 



6S1/2 


-190.88 


-128.28 


6S1/2 


-51.16 


-35.99 


6Si/ 2 


-169.96 


-123.44 


6Si/ 2 


-468.88 


-369.62 


6Si/ 2 


-90.62 


-73.46 




6p 



6p 



7s 



7P,. 



FIG. 2: The ratios of pair-correlation and core-polarization 
effects w.r.t. the DF values 



TABLE IV: Contributions of the 6s!/2 core electron (in MHz) 
to the core-polarization effect using the RMBPT(2) approxi- 
mation (third column) and RCC theory (forth column) 



states (which were absent in the DF approximation). 
It is interesting to note that the DF values for these 
two states deviate from their respective experimental 
values in opposite directions, so that the sign of the 
correlation contributions are opposite for the two cases. 
This is further supported by the results based on 
second order relativistic many-body perturbation theory 
(RMBPT(2)) which are given in table III. Here the 



Terms 6p 1/2 6p 3/2 7s 1/2 7pi/ 2 

state state state state 



O 11513.5 918.3 7822.9 1983.1 

O-O 665.3 -43.7 983.3 85.4 

OS lv + cc 952.2 78.4 2122.6 326.6 

OS 2v + cc 1188.2 -591.0 1916.8 35.6 

S\ v dSi v 21.0 1.6 164.6 14.1 



S\ v OS 2v + cc 22.2 0.6 180.2 19.2 
Sl v OS 2v + cc 149.6 194.61 298.8 18.7 

Important effective two-body terms of O 



Sl v OT! + cc -20.2 2.0 14.6 -0.76 
S f 2v OT 2 + cc -160.2 -12.6 -135.4 -21.64 
Norm. -88.5 -6.7 -181.8 -22.98 

Total 12903.7 623.2 11158.6 2263.5 
Experiment 13000 583(21) 



TABLE V: Contribtions of different coupled-cluster terms to 
the Pb + hyperfine stucture constant, cc stands for the com- 
plex conjugate part of the corresponding terms 



dominant contributions to RMBPT(2) as shown in Fig. 
B are given explicitly. While electron correlation at this 
level is substantial for all the states, it is dramatic in 
the case of 6p 3 / 2 state because of the unusually large 
and negative core-polarization. In table IV we give the 
significant contributions to the core polarization which 
arises from the interaction of the outermost core 6s and 
the valence 6^3/2 electrons (third column). The sum of 
these individual contributions is -971.5 MHz and after 
taking into account the polarization of all the other core 
electrons, a net contribution of -814.6 MHz is obtained. 
The tremendous size of this second-order correction 
suggests that an all order method like coupled-cluster 
theory is necessary for a correct quantitative description 
of the correlation effects in Pb + . This is indeed reflected 
in the results given in table V. Again, the 'bare' operator 
O represents the DF approximation, i.e. excluding any 
correlation effects. Several important correlation con- 
tributions for our RCC calculations are also presented 
in table V. O = e T+ Oe T denotes the 'dressed' operator 
containing the contributions of the closed-shell cluster 
amplitudes in Eq. i|U.6|) . Although core-polarization 
(OS2V) and pair correlation (OSiv) are the dominant 
correlation effects, core correlation effects (O — O) are by 
no means negligible; they amount to about 9% for the 7s 
state. Summing up all the contributions given in table 
V leads to significant improvements in our calculated 
values of the hyperfine constants of the 6pi/2 and 6p 3 / 2 
states. The sub one percent (0.7%) agreement between 
the former and experiment is indeed spectacular. A 
similar agreement cannot be expected for the latter state 
which is characterised by extremely peculiar correlation 
effects. Even so, the hyperfine constant for this state 
differs from experiment (3.6% error bar) by a little less 
than 7%. This is certainly remarkable considering that 
the corresponding discrepancy at the level of RMBPT(2) 
is 48%. It is interesting to note from table IV (fourth 
column) that the core-polarization contributions in RCC 
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theory follow the same trend as in RMBPT(2). 

The plot in fig. 2 highlights the relative importance 
of the core-polarization and pair correlation for the 
different states. It is instructive to point out that, unlike 
the hyperfine constant in the ground state of Ba + @, 
core-polarization effects are larger than pair-correlation 
for the ground and first excited state, i.e. 6pi/2 and 

6p ; 



3/2 



states of Pb^ 



This is the result of the much 



stronger valence-core interactions in Pb + compared to 
Ba+. 

In summary, the strength of RCC theory has been 
exploited to obtain for the first time a quantitative 
understanding of the interplay of relativistic and 
correlation effects in the magnetic dipole hyperfine 
constants for Pb + . It has been demonstrated that 
the results of the DF and RMBPT(2) approximations 
differ substantially from the measured values of the 



hyperfine constants. However, the inclusion of single, 
double and a subset of triple excitations to all orders 
in the framework of RCC theory leads to a dramatic 
improvement in the results. The relevance of the present 
work extends beyond hyperfine interactions in Pb + . 
Our results highlight the fact that a judicious use of 
RCC theory can yield accurate results for properties 
that are sensitive to the nuclear region. Indeed, this has 
important implications for Tl, which like Pb + is a heavy 
atomic system with strongly interacing configurations 
and is one of the leading candidates for the study of 
parity nonconservation due to neutral weak currents 

We are grateful to Prof. Werth for his valuable discus- 
sions and suggestion for this calculation. It was possible 
to contact with Prof. Werth through DST-DAAD ex- 
change programme. The calculation is carried out using 
the Tera-flopp Supercomputer in C-DAC, Bangalore. 
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